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Abstract 

We perform a detailed study of the importance of loop corrections when deriving bounds on 
heavy-active neutrino mixing in the context of general Seesaw mechanisms with extra heavy right- 
handed neutrinos. We find that, for low-scale Seesaws with an approximate B — L symmetry 
characterized by electroweak scale Majorana masses and large Yukawas, loop corrections could 
indeed become relevant in a small part of the parameter space. Previous results in the literature 
showed that a partial cancellation between these important loop corrections and the tree level 
contributions could relax some constraints and lead to qualitatively different results upon their 
inclusion. However, we find that this cancellation can only take place in presence of large violations 
of the B — L symmetry, that lead to unacceptably large contributions to the light neutrino masses 
at loop level. Thus, when we restrict our analysis of the key observables to an approximate B — L 
symmetry so as to recover the correct values for neutrino masses, we always find loop corrections 
to be negligible in the regions of the parameter space preferred by data. 
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I. INTRODUCTION 


The origin of the observed pattern of neutrino masses and mixings in neutrino oscillation 
experiments (see e.g. Ref. [1] for a recent summary) comprises one of the few experimental 
evidences for physics beyond the Standard Model (SM) of particle physics. The simplest 
and most popular extension to account for these experimental observations consists in the 
addition of right-handed neutrinos to the SM particle content. Given their singlet nature, 
a Majorana mass term for the right-handed neutrinos is directly allowed in the Lagrangian, 
thus inducing a new mass scale -the only one unrelated to electroweak (EW) symmetry 
breaking- to be determined by data. Depending on the size of this scale its phenomenological 
consequences are very different. 

One of the most appealing choices is that this new Majorana scale is high, leading to the 
well-known Seesaw mechanism [2-5] and providing a rationale for the extreme smallness of 
neutrino masses when compared to the rest of the SM fermions and the EW scale. Values 
for the neutrino Yukawa couplings ranging between that of the electron and that of the 
top quark would lead to Majorana masses between the EW and the grand unification scale. 
Unfortunately, even for the lightest mass choice, any phenomenological consequence beyond 
neutrino masses tends to be hopelessly suppressed if the extra degrees of freedom only couple 
to the SM through their Yukawa interactions. 

However, the smallness of neutrino masses could derive from symmetry arguments [6-9] 
rather than a hierarchy of scales. Indeed, the Weinberg operator [10] leading to neutrino 
masses in Seesaw mechanisms is protected by the B — L symmetry, conserved in the SM 
and violated in two units by the Weinberg operator. Thus, if the pattern of the Yukawa 
couplings and Majorana masses in a Seesaw realization is such that it conserves B — L, 
the Weinberg operator will never be generated and the SM neutrinos will remain massless, 
even for Y v ~ 1 and Majorana masses of the order of the EW scale. Small violations of 
B — L in this pattern would thus induce the small neutrino masses observed in oscillation 
experiments. In this class of models fall the popular inverse [6, 7] or linear [11] Seesaw 
mechanisms which, contrary to the canonical type-I Seesaw, would lead to an extremely rich 
phenomenology through the large mixing allowed between the new extra sterile neutrinos 
and their SM siblings implying observable contributions in lepton flavour violating (LFV) 
processes, universality violation and signals in electroweak precision observables. It is then 
of interest to fit all these available data to determine the allowed values of the mixing of the 
heavy neutrinos with the SM charged leptons, examples of these constraints can be found 
in Refs. [12-29]. 

When deriving such constraints on heavy-active neutrino mixing, it was recently pointed 
out in [24] that loop corrections involving the extra heavy neutrinos played an important 
role, obtaining qualitatively different results to those derived by staying at leading order. In 
particular, it was shown that corrections to the T parameter [30, 31] could be sizable and 
that these, in turn, would affect the determination of Gp through /i decay competing with 
the tree level effects. Since the value of Gp from p decay is generally in good agreement 
with the measured value of Mw and other determinations of sin dye, in [24] it was found 
that the constraints stemming from these datasets could be weakened at loop level through 
partial cancellations between the tree level corrections and the T parameter contribution. 
Furthermore, the invisible width of the Z, which is in slight tension with the SM prediction, 
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is modified at tree level through the presence of extra heavy neutrinos, while the oblique 
corrections computed in [24] were found to be subleading. Thus, by accounting for these loop 
corrections, good fits with relatively large heavy-active mixing were found in [24], since it is 
possible to alleviate the tension in the invisible width of the Z without seriously affecting 
the determination of Gp in /i decay through the partial cancellation of the tree and loop 
level contributions. 

However, when Ref. [25] also investigated the relevance of the T parameter the same can¬ 
cellation was not reproduced and in [27] it was argued that loop contributions should always 
be negligible, since the heavy-active mixing that controls the strength of the couplings of 
the new degrees of freedom is bounded to be small ( 6 2 < 10~ 2 ). Therefore, new tree-level 
bounds were derived instead through more updated fits to available data. While this argu¬ 
ment is generally true, models based on an approximate B — L symmetry are characterized 
by large Yukawas and EW-scale Majorana masses, thus, even if loop corrections through 
weak interactions further suppressed by 6 2 are indeed negligible for all practical purposes, 
when the loop corrections are mediated by heavy neutrinos and/or the Higgs held or its 
Goldstones, the coupling involved in the vertex is no other than the large Yukawa coupling, 
so that loop corrections can indeed become relevant, as stated in [24], However, not only 
the oblique corrections computed in [24] fall in this category, since the effect of the large 
Yukawa interactions does not vanish in the limit of massless neutrinos and charged leptons. 
Indeed, some vertex and box corrections involving large Yukawas are found not to vanish in 
the massless limit for light leptons (see e.g. [32]). 

In this work we clarify the importance of loop contributions to the determination of the 
heavy-active neutrino mixing including all loop corrections mediated by the potentially large 
Yukawa interactions. We find that, as discussed by [24], the most relevant of these corrections 
are indeed the ones encoded through the oblique parameters but, in order to make them 
competitive with the tree-level contributions, EW scale Majorana masses and Yukawas on 
the very border of perturbativity are simultaneously required. Furthermore, we find that, as 
long as B — L is conserved, the T parameter is always positive, so that the partial cancellation 
discussed in [24] cannot take place in such a setup. Large violations of B—L are thus required 
to achieve the negative and sizable values of T capable of reproducing the effect. But these 
large violations of B — L would render the Weinberg operator unprotected and, in presence 
of the EW-scale Majorana masses and large Yukawas required for T, radiative corrections 
lead to unacceptably large contributions to the light neutrino masses, even if these are fixed 
to their correct value at tree level by means of the Casas-Ibarra parametrization. Thus, 
when we impose an approximate B — L symmetry with only 3 extra heavy right-handed 
neutrinos, we always find that loop corrections are irrelevant when deriving bounds on the 
heavy-active neutrino mixings. 

This paper is organized as follows: In Section II we detail the parametrization employed 
for our study. In Section III we list the observables we analyze in our global fits. In Section IV 
we present our findings and discuss the importance of loop effects in the global fits as well 
as the necessity of large violations of B — L in order to obtain partial cancellations between 
the tree and loop level contributions. Finally, in Section IV we summarize our results and 
present our conclusions. 
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II. PARAMETRIZATION 


In this work we explore the constraints that can be derived through various EW ob¬ 
servables on the extra neutrino mass eigenstates mixing with charged leptons in a Seesaw 
scenario: 

& = - -N&iMNhNg - (Y N ) ia Ni<t>^Z + H.c. . (1) 


Here, (j) denotes the SM Higgs field, which breaks the EW symmetry after acquiring its 
vev uew- We have also introduced the Majorana mass Mn allowed for the right-handed 
neutrinos N ^ as well as the Yukawa couplings between the neutrinos and the Higgs field. 
We will restrict our study to the extension of the SM by 3 right-handed neutrino fields. The 
vev of the Higgs will induce Dirac masses mp = uew^ 7 /v/v / 2. Thus, the full 6x6 mixing 
matrix U is the unitary matrix that diagonalizes the extended neutrino mass matrix: 


/ 0 m'jj 
\m D M n 


U 


( m 0 
0 M 


( 2 ) 


where m and M are diagonal matrices containing respectively the masses of the 3 light ty 
and 3 heavy N t mass eigenstates. The diagonalizing matrix U can be written as [33]: 


where 




/ 
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( 3 ) 


( 4 ) 


and Dpmns is, approximately , the PMNS matrix measured in neutrino oscillation experiments 
up to the non-Unitary (Hermitian) corrections from c. For alternative parametrizations of 
the full mixing matrix see Refs [34-38]. Indeed, due to this Hermitian correction, the 
actual PMNS matrix appearing in charge current interactions mixing the light neutrinos 
and charged leptons will, in general, not be Unitary and we will refer to it as N: 


N — cUpmns (5) 

The general matrix 0, representing the mixing between active (z/ e , z/ M and u T ) and heavy 
(W, N -2 and W) neutrino states, and the mass eigenstates m and M are determined from 
Eq. (2) which leads to: 

C ^PMNS m ^PMNS C = ~ s Ms*. (6) 

In the Seesaw limit, that is M^ mp, these conditions reduce to the well-known results: 

0 ~ m^M^ 1 

^PMNS m ^PMNS — ~ mt D^N lm D = 

M ~ M n . (7) 
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Notice that, naively, the mixing between the active and heavy neutrinos 00l ~ m/M 
and, given the smallness of neutrino masses m, the mixing effects we will study here 
would be unobservably small. However, in the context of Seesaw mechanisms with an 
approximate conservation of B — L such as the inverse [6, 7] or the linear [11] Seesaws, 
this symmetry suppresses the neutrino mass m while allowing a sizable mixing. This ap¬ 
proximate symmetry not only ensures an equally approximate cancellation in the combina¬ 
tion rrijjM^rriD leading to the observed neutrino masses while allowing large -potentially 
observable- 0©t = but also ensures the radiative stability and technical natu¬ 

ralness of the scheme [39]. 

When extending the SM Lagrangian by only 3 new singlet (right-handed neutrino) fields 
essentially the only neutrino mass matrix with an underlying L symmetry that leads to 3 
heavy massive neutrinos is [40] (see also Ref. [41]): 


m D 


vew 

75 


f Y Y Y \ 

1 e 1 fi 1 t 

eiYJ ei Y' 

V ^ Y " e 2 Y” e 2 y" 


and M n = 


1 l-i \ A /x 3 ^ 

A /1 2 /l4 1 

^ h3 h4 A' j 


( 8 ) 


with all and fij small lepton number violating parameters (see also Ref. [42] for a particular 
scenario where these small parameters arise naturally). Indeed, setting all e, = 0 and fi 3 = 0, 
lepton number symmetry is recovered with the following L assignments L e = = L T = 

Li = — L 2 = 1 and = 0. In Eq. (7) this leads to: rh = 0 (3 massless neutrinos in the 
L-conserving limit), M\ = M 2 = A (a heavy Dirac pair) and M 3 = A' (a heavy decoupled 
Majorana singlet), but: 
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^ew 

~2A 


( -%Y* Y* 0^ 
-iy; Y* 0 
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— id e 9 e 0 

—i0,j 9fj, 0 
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and 00^ = 


( \9e\ 2 9 e o; 6 e 9* ^ 
A 1/3 |2 n A* 

v m; m; %\ 2 j 


(9) 


Thus, vanishing light neutrino masses can still be associated with arbitrarily large mixing 
between the heavy Dirac pair and active neutrinos and, for these kind of Seesaw scenarios, the 
bounds on the mixing we will explore are complementary and independent to the stringent 
constraints on the absolute light neutrino mass scale. 

The small L- violating parameters e* and /i ? will induce small non-zero neutrino masses and 
mixing among these light mass eigenstates but will only translate in negligible perturbations 
to the matrix 0. With the simple form in Eq. (9) for the heavy-active mixing, the series 
expansions in Eq. (4) can be added exactly obtaining: 


sin 6 

s = —— 0 


and 


c — I — 


1 — cos 9 
02 


00 f , 


( 10 ) 


with 

9 = \J\9 e V + |$^| 2 + |$t | 2 ■ (11) 

Regarding the role of the e % and fij parameters in the generation of the light neutrino 
masses and mixings observed in neutrino oscillations, all of them except /ii and /i 3 will lead 
tom^O through Eq. (7) when switched on: 
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e 2 — ( m D A lm D + m D A' 1 m D ) , 


( 12 ) 


with 


m D = 


^ew 

7 f 


(X„Y„Y r 


Ill 


D — 


^EW 

V 5 


(Y:,Y'K) and 
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D — 


^EW 


wxx)- 


i (13) 

Indeed, even though /i 4 and /13 do violate L, upon their inclusion the mass matrix in Eq. (2) 
does not increase its rank, which, in absence of the other e, and fi 3 , is only 3 and thus 3 
massless eigenstates are still recovered 1 . The parameters /i 2 and yU 4 do contribute at tree 
level to generate light neutrino masses, however, their effect can be absorbed in a redefinition 
of the vectors rri^ and mj, as follows: 


/ / AD 

eim D eim D - —m D 


and 


e 2 mD e 2 mj) 


/i-4 

A m ° 


(14) 


up to contributions with two extra powers of the small L-violating parameters. Thus, in 
presence of non-zero it is enough to consider their contribution to the generation of 
neutrino masses which reads: 


m = eim^A 1 mD + eim^A 1 m / D + e^m^A' 1 m / ] ( ) . (15) 

Notice that the last term in Eq. (15) is suppressed by two powers of e 2 while the others only 
by one power of e 4 . However, e 2 (and /13 and /i 4 ) violates L by one unit while e 4 (and fi\ 
and /r 2 ) by 2. Hence, if the source of L-violation is by one unit it is expected that e 4 ~ e\. 
Thus, for full generality, we will keep the last term in Eq. (15). The six free parameters 
encoded in and allow to give mass to the three mass eigenstates observed in neutrino 
oscillations as well as the possibility of reproducing any mixing pattern including the, yet 
unknown, CP-violating phases of Dirac and Majorana types encoded in the PMNS matrix, 
while leaving m^, and hence 0, s and c, mostly unconstrained 2 . One of the three elements 
of mo is, however, fixed by the other two, the values of the light mass eigenstates and the 
elements of the PMNS matrix when solving for Eq. (15) obtaining the following relation: 

Yr - —y-——v— (Y e (rhe^rh^ - fh eT fhpn) + 

mj M - rn ee m Nl 

Y„ ( m efl m eT - rheeTh^r) - ^Y^rh^ - 2Y e Y ll m eil + Y^m ee x ( 16 ) 

x yjm\ T m^ - 2m e/1 m eT rh flT + rh ee m + m 2 efl m TT - rh ee rh^rh T ^j , 


1 Notice that, even if yi and ^3 do not induce neutrino masses at tree level, the L symmetry protecting 
them is now broken and loop contributions would appear instead [43]. 

2 In contrast, neglecting the last term in Eq. (15) would lead to the more constrained scenario explored 
in detail in Ref. [44], with a massless neutrino and a mixing pattern in 0, s and c determined up to an 
overall factor from the observed neutrino oscillation parameters. This scenario has also been studied in 
Refs. [45-49] 
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Parameter 

\Y e \ x |Y m | 

\Y e \ -\Y „| 

m\ [eV] 

A [GeV] 

Phases: a e , a^, 6, ctu & «2 

Osc. data 

Range 

( 0 ,IQ" 4 ) 

(- 0 . 1 , 0 .!) 

(io- 5 ,1) 

( 10 3 , 10 4 ) 

(0, 2 tt) 

fixed [ 1 ] 


TABLE I: The 9 free parameters of our scan: the modulus and phase of the electron and muon 
Yukawas \Y e \, a e and a the Majorana mass scale A, the absolute neutrino mass mi and 
the 3 yet unknown CP-violation phases (Dirac and Majorana) in the PMNS mixing matrix: 5, a± 
and a. 2 - The PMNS mixing angles and mass splittings are fixed to their best fit from the global 
analysis in Ref. [1]. 


where rh = — ^p MN s m ^p MN s * s the mass matrix of the flavour eigenstates. Thus, in our 
numerical exploration of the parameter space in Section IV we will consider the 9 free 
parameters summarized in Table I. 

An alternative parametrization extensively used in the literature is the so-called 
Casas-Ibarra parametrization [50]. This parametrization introduces the matrix R = 
iM~ 1 ^ 2 m D UpM'Ns rn ^ 1 ^ 2 exploiting the fact that, from Eq. (7), R has to be (complex) or¬ 
thogonal. The main advantage of this parametrization is the ability to easily recover the 
Yukawa couplings through the heavy mass eigenvalues M and the low energy observables 
b^pMNS an d m together with the elements of R as rri£, = —fM 1 // 2 i?m 1 // 2 f/p MNS . However, 
the physical range of the parameters contained in R can be cumbersome and a physical 
interpretation of their values is not immediately transparent, see Ref. [51] for a detailed 
discussion. Moreover, these relations only hold at tree level 3 . Thus, when values of R are 
chosen so as to allow sizable low energy phenomenology through large Yukawas and low M, 
it is important to check if the pattern displays an approximate B — L symmetry. Otherwise, 
loop corrections to the unprotected Weinberg operator, that is to Hpmns and m, will exceed 
present constraints even if their values were correct at tree level. For this reason we rather 
chose to perform the scan through the parameters summarized in Table I. 

At energies much below the masses of the heavy neutrinos A and A' the effects of their 
mixing O manifest dominantly through deviations from unitarity of the lepton mixing matrix 
N. Since any general matrix can be parametrized as the product of an Hermitian and a 
Unitary matrix, these deviations from unitarity have been often parametrized as [53]: 


N — (I — 77)UpMNS (17) 

where the small Hermitian matrix r] (also called e in other works) corresponds to the coeffi¬ 
cient of the only dimension 6 operator obtained at tree level upon integrating out the heavy 
right-handed neutrinos in a Seesaw scenario [54] and, in our parametrization it would be 
given from Eqs. (5) and (10) by: 

, = izX^ee. (18) 


3 See Ref. [52] for a generalization of the Casas-Ibarra approach to loop level. 
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III. OBSERVABLES 


In this section we introduce the list of observables used for our analysis. While a more 
comprehensive set could be considered (see for example Ref. [27]), we have rather chosen 
the most representative of these observables since extending the analysis to the loop level 
for the whole set would be cumbersome and the dominant constraints as well as the main 
effects pointed out in [24] are contained in a smaller subset. We will thus present both the 
1-loop contributions and the experimental constraints for a total of 13 observables. The loop 
amplitudes of the processes have been computed exploiting the Goldstone-boson equivalence 
theorem [55] under the assumption that the mass of the extra neutrinos M % is larger than the 
gauge boson masses; i.e. Mi > Mw,z■ Thus, we have made the simplifying assumption that 
the most relevant loop corrections are those were the loops are mediated by either the Higgs 
boson, h , the Goldstone bosons and </>° or the heavy Majorana neutrinos. Indeed, this 
forces the vertexes to involve the potentially large Yukawa couplings (the only couplings 
that can be relevant at the loop level) and the corrections from including the transverse 
components are suppressed by M 2 vz /M 2 r . The set of 13 independent observables analyzed 
in this study is composed of: 

• 8 ratios constraining electroweak universality: RJ iel R ^ fl , R? v e , R}^, R^ e , R 1 ^, 7?J /e , R} 

• The invisible Z width 

• The W mass M w 

• 3 rare flavour-changing decays: /i — > ey, r —> /ry and r —> ey 

All of them will be determined as a function of the three most precise electroweak measure¬ 
ments: a, Mz and G M ( Gf as measured from /i decay) [56]: 

a = (7.2973525698 ± 0.0000000024) x 10" 3 , 

M z = (91.1876 ±0.0021) GeV, (19) 

= (1.1663787 ±0.0000006) x 10 ~ 5 GeV" 2 . 

All observables will receive contributions from the loop corrections to the W and Z boson 
propagators through the diagrams in Fig. 1 . These contributions are encoded in the flavour- 
universal corrections Sy^z th a t can be found in Eq. (67) in the Appendix. We now list the 
further corrections exclusive to each of the observables considered: 


A. /x decay, Gf and My y 


Our input value for Gp is determined through /i decay, but this process will receive 
corrections both at the tree and the loop level (see Fig. 2). Thus, the value determined from 
/i decay, G is related to Gp by: 


= mlG 2 F 

M 192vr 3 


(l-|0 e | 2 -|^| 2 ±2^ ivJV ±5G) = 


rrjCtl 

192tt 3 ’ 


( 20 ) 






FIG. 1: 1-loop correction of the new heavy neutrinos to W and Z propagators. 



FIG. 2: 1-loop corrections to /i decay. 

with 

5G = 2 Re[VY + + 5 e CT w + <fp T w * + B Me ] (21) 

and where N is the flavour-universal W propagator correction, hp T w and V/' are the 
flavour-dependent lepton propagator and vertex contributions (see Eqs. (55) and (57) in 
the Appendix), and £> Me encodes the box diagram contribution computed in Eq. (59) in the 
Appendix. 

From Eq. (20), we find: 

Gl = G 2 f ( 1 - \d e \ 2 - \9,\ 2 + 25^ ivN + 5G) . (22) 

The second and third terms in Eq. (22) correspond to the tree level correction, the 
fourth term is the universal 1-loop oblique correction which is given in Eq. (67) of the 
Appendix. This particular expression, when used in an observable mediated by the Z and 
thus corrected through 2 8^ mv N , leads to a common correction to these observables given by 
1 — \9 e \ 2 — I^J 2 — 2 aT (see Eqs. (67) and (63)). This common dependence on the tree level 
and oblique corrections is the source of the cancellation analyzed in Ref. [24], 


9 










The the W mass is also correlated to G F through 



tv a 

s/2G f s*,(1- Ar)’ 


( 23 ) 


with Ar = 0.03639 =F 0.00036 ± 0.00011 [56]. Thus, the corrections induced at both the tree 
and loop levels by the heavy neutrinos from Eq. (22) can be probed by the measurement of 
M w in LEP and Tevatron [56]: 


M w = 80.385 ± 0.015 GeV. 


(24) 


B. Invisible Z width 


The determination of the number of light active neutrinos by LEP through the invisible 
width of the Z provides a constraint to heavy neutrino mixing already at the tree level. 
Additional loop corrections are induced through the diagrams in Fig. 3 which lead to: 


r. 

L inv 


3 


E 


GpMzP 

24^277 


(Zij + Zji ) , 


(25) 


where p encodes the SM loop corrections to the process and 


Zii = |Cbf(l + ^ niv ) + 2 Re[C% (5§ t z + V?) ] , 


(26) 


with 


Ca = 


E u*.u ■ 

u m u O‘3 


a=e,ii,T 


(27) 


and S^ T z and V z the lepton and vertex corrections shown in Eqs. (56) and (58) in the 
Appendix. 

Eq. (25) is often used to determine the number of active neutrinos N v lighter than Mzj 2 


as: 


_ G F M\pN v 
Y 2\[2 


The measurement by LEP of r inv = (0.4990 ± 0.0015) GeV combined with Eq. (28) leads 
to [56]: 

N v = 2.990 ±0.007. (29) 


We will exploit this result together with Eq. (25) to derive constraints on CV, and, hence, 
on the heavy neutrino mixings. 


C. Universality ratios 

Electroweak coupling universality is strongly constrained through ratios of lcptonic decays 
of K, 7r, W or charged leptons. In these ratios many uncertainties cancel and a clean 
constraint can be derived. These observables are corrected both at the tree and loop level, 
for instance, BZ = Y (tt~ —$■ /Y ( ti~ —$■ eV e ) is corrected by the diagrams in Fig. 4. 
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BR (7r + — > e + v e ) 

(1.230 ±0.004) x 10 -4 

BR(n + 

(99.98770 ± 0.00004) % 

BR ( T~ TT-fr) 

(10.83 ± 0.06) % 

BR (K + e + v e ) 

(1.581 ±0.008) x 10“ 5 

BR{K+^ p+vj 

(63.55 ± 0.11) %10 -5 

BR(t~ -> K~u t ) 

(7.00 ±0.10) x 10" 3 

BR(W + -> e + v e ) 

(10.71 ±0.16)% 

BR(W+^^ + ^) 

(10.63 ±0.15)% 

BR (W + -> r+i/ T ) 

(11.38 ±0.21)% 

BR(t~ —>• 

(17.41 ± 0.04) % 

Bi? (r“ —> e~V e u T ) 

(17.83 ± 0.04) % 

^7T± 

(2.6033 ± 0.0005) x 10" 8 s 


(1.2380 ±0.0021) x 10“ 8 s 

T t 

(290.3 ±5.0) x 10“ 15 s 


(2.1969811 ±0.0000022) x 10" 6 s 

m n ± 

139.57018 ± 0.00035 MeV 

m K ± 

493.677 ± 0.016 MeV 

M\y 

80.385 ± 0.0015 MeV 

m e 

0.510998928 ±0.000000011 MeV 

m fj. 

105.6583715 ± 0.0000035 MeV 

m T 

1776.82 ±0.16 MeV 

SR^e 

(-0.374 ±0.001) 

SR^r 

(0.0016 ± 0.0014) 

sr%t 

(0.0090 ± 0.0022) 


TABLE II: Input values used for the constraints on weak universality from ratios of meson and 
charged lepton decays. 
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2 


r — 

1 inv — 



FIG. 3: 1-loop corrections to the invisible decay of the Z. 



Thus, the general expression for the ratio of lepton flavours a and (3 is given by: 


„ p SM l-|W 2 + 2Re[vy + e"'] 
** aS 1 - 1 0f, f 2 + 2Re [ vy + Sf w ] ’ 
where is the SM value for this ratio, for example, for tt decay: 

r ,sm = ( ( m l - m l) 1 1 

\">J (ml -ml)) 1 + snip 


(30) 


(31) 


and where SR*g are the SM radiative corrections to this process [57]. Notice that the 
flavour-universal contributions from the W propagator cancel in the ratio. 

The predicted values of these ratios are computed through Eqs. (30) and (31) with data 
form [56, 58] and compared to the experimental measurements of the decay rates in our 
global fit. This data is summarized in Table II. 


D. Rare decays 

The presence of extra heavy neutrinos beyond the three light ones participating in low 
energy weak processes induces deviations from unitarity in the PMNS matrix. Thus, the 
GIM cancellation [59] suppressing flavour-changing processes does not take place and strong 
constraints on the presence of these extra neutrinos can be derived. Moreover, the extra 
heavy neutrinos themselves also mediate the flavour-changing processes, such as radiative 
leptons decays l a —>■ in Fig. 5. The contribution from both the heavy and light neutrinos 
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FIG. 5: Extra neutrino contributions to the fi —y e 7 decay. 


is given by: 


T (Z a —?• Z^y) 3a 


UotkU\pF(xk) 


k =1 


r (Z Q -a Z^ Q ^) 32vr {UW) aa (UW) 


pp 


Ml 


where x k = jpS and F(x k ) is given by: 


^Ofc) 


1 i) — 43x fc + 78xl — 49xl + 4xi + 18x? In x k 


3(x k - l ) 4 

Thus, for heavy neutrino masses much larger than M\y\ 

r (Z Q —y Ip y) 3a 


\e a e;\\F(oo)-F(o)y. 


(32) 


(33) 


T (l a —y IpVaVp) 32n 
The prediction from Eq. (32) will be compared with the existing upper bounds from [56]: 


(34) 


BR, lie < 5.7 x 10~ 13 , 
BR re < 3.3 x 10 ~ 8 , 
BR, ril < 4.4 x 1CT 8 . 


(35) 

(36) 

(37) 


Notice that these bounds are quoted at the 90% CL so they will be rescaled to la to build 
the corresponding contribution to the y 2 function. 


IV. RESULTS 

A. Constraints from the global fit 

With the 13 observables discussed in Section III we build a x 2 function depending on the 9 
parameters listed in Table I. Given the large dimensionality of the parameter space, we make 
use of Markov chain Monte Carlo (MCMC) techniques for efficient parameter exploration. 
In particular, we implement importance sampling based on the Likelihood obtained from 
the observables through a Metropolis-Hastings algorithm. The range in which the 9 free 
parameters are varied is also summarized in Table I. We have run simultaneously 5 different 
chains through the MCMC algorithm and have verified that good convergence (better than 
R — 1 < 0.035 [60]) for all parameters has been achieved. The results of the runs thus 
provide a good sample of the y 2 values in the preferred regions of the parameter space and 
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FIG. 6: Contours for 9 e , 6 M and 0 T at la (red), 90% CL (black) and 2 a (blue). The left panels are 
obtained for normal hierarchy and the right for inverted. 


have been used to marginalize over different subsets of the model parameters. In this way, 
we will present 2D and ID frequentist contours on the more phenomenologically relevant 
parameters of the model. The post-processing of the chains to derive the allowed confidence 
regions has been performed with the MonteCUBES [61] user interface. 

In Fig. 6 we show the results of our MCMC scan for the 2 degrees of freedom constraints 
of different combinations of the heavy-active mixings 9 a defined in Eq. (9). The contours 
correspond to the la, 90% and 2a frequentist confidence regions. The upper panels show the 
bounds in the two combinations we choose to more directly sample (see Table I): \9 e \ x \9^\ 
and 1 9 e | — l^j. The rationale behind this is apparent upon inspection of Fig. 6 . Indeed, 
the constraints on the product are more than one order of magnitude smaller than those 
derived from the difference of the couplings \/\9 e \ x \9^\ -C \\9 e \ — \9 fl \ \, leading to a very 
pronounced hyperbolic degeneracy in the panels of the middle row, which contain the same 
information directly depicted as a function of 9 e and 9 )1 . Thus, this particular choice of 
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sampling parameters allowed to scan the hyperbolic degeneracy much more efficiently and 
speed the convergence of the MCMC. This very strong constraint in \9 e \ x 1 9 t ,| stems from 
the strong bound on jjl —> ey from MEG that, from Eq. (34), sets a very stringent limit on 

M- 

Finally, the lower panels of Fig. 6 contain the constraints derived for the mixing with 
the r flavour 9 r . Notice that Y T , and hence 0 T , was not a free parameter of the fit but was 
rather obtained from the other two Yukawas and the light neutrino masses and mixings from 
Eq. (16). This is the source of the observed correlation between the values of 9 e and 9 r . 
Notice also that, since the particular pattern of light neutrino masses plays an important 
role in Eq. (16), the left (normal hierarchy) and right (inverted hierarchy) panels of Fig. 6 
display different correlations. 

In Fig. 7 we show the individual constraints that can be derived on 9 e , 6% and 9 r (from 
top to bottom) for a normal (left) and an inverted (right) hierarchy after marginalizing over 
all other parameters. We generally fold a slight improvement of the fit to the observables 
considered when some amount of mixing is present. In particular, we fold that non-zero 
mixing with the electron is preferred at around the 90% CL by our dataset. Mixing with the 
tau flavour is also favoured for normal hierarchy due the correlations implied by Eq. (16). 
At the lcr level, mixing with the /i flavour is significantly constrained due to the preference 
of some universality bounds (from 7r and r decays) for a slightly reduced coupling to the 
electron with respect to the muon. Thus, since universality constraints are corrected by 
1 — \9 a \ for each flavour, a non-zero 9 e is preferred in the fit while 9 M is kept at small values 
to satisfy the constraint from /i —>• ey. Beyond the la level, the mixing with the electron 
is allowed to become small and thus the constraint on /i mixing at 2cr is much weaker than 
naively expected from the la region. The limits of the 1 and 2cr regions for the three mixing 
parameters are summarized in Table Ilf. 



Be 


Or 


la 

2 a 

la 

2 a 

la 

2(7 

NH 

D 034+ 0 - 009 

u.uj4„ 0 .o 14 

< 0.050 

< 3.2 • 10" 4 

< 0.037 

0.018±g;gii 

< 0.049 

IH 

D 035 +0, ° 09 
U.UoO_ 0 0 44 

< 0.051 

< 3.3 • 10" 4 

< 0.037 

< 0.031 

< 0.044 


TABLE III: Constraints on 9 e , 8 and 0 T for normal and inverted hierarchy. 

In Fig. 8 we show a comparison of the breakdown of the contributions of the different 
observables to the total y 2 f° r the SM (left panel) and our best fit (middle panel) as well as 
the difference of the two (right panel). It can be seen that some of the existing tension of 
the SM with the invisible width of the Z can be alleviated by the presence of heavy neutrino 
mixings and also the agreement between the kinematic determination of M\y and its SM 
value from Gf, a and Mz is improved. As already discussed, the universality constraints 
from 7r and r decays are also in better agreement when some mixing with the electron is 
present. On the other hand, universality tests from kaon decays rather point in the opposite 
direction. Thus, at the end, the preference for non-vanishing heavy-active mixing is mild 
and the final improvement of the y 2 with respect to the SM value is 3.7, not quite reaching 
the 2 a level. Notice that, even if the number of free parameters in the fit is rather high, the 
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FIG. 7: Ay 2 (marginalized over all other parameters) for 9 e , 6^ and 6 r . Left panels show results 
for a normal hierarchy and right panels for inverted hierarchy. 

observables actually depend on the combinations \9 e \, \ 0 /t | and \9 r \ only (and A when loop 
corrections are relevant). Thus, the reduction by 3.7 of the y 2 should be attributed to the 
introduction of 3 (or 4) new parameters rather than 9. 

Regarding the importance of the loop effects considered, we have performed a second 
set of MCMC runs where all loop corrections have been removed. The results of these 
simulations are essentially identical to the ones stemming from the full computation. By 
adding to the chain output also the value that the T parameter took in the simulations, 
we find that its preferred values are ~ 10“' — 10” 6 , negligible with respect to the best fit 
values of the tree level contributions. In order to understand this apparent lack of relevance 
of the loop corrections and the T parameter in particular, in direct contrast to the results 
presented in [24], we will now analyze in further detail the regions of the parameter space 
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FIG. 8: Contributions from the different observables to the y 2 . Left plot shows the SM values. 
Middle plot shows the contributions from three right-handed neutrinos in the best-fit point. Right 
plot shows A x 2 = Xi(SM) — Xi(BF) for every observable i. 

in which T could be relevant and the necessary conditions for the cancellation with the tree 
level contributions to take place. 


B. The T parameter 


The leading contributions (not suppressed by the light neutrino or charged lepton masses) 
to the T parameter are given by [24]: 


aT = 


a 




(e K 


UajUpUhfr + u:,u aj u;,Lx gij ) 


(38) 


where 

, M?MJ M, , 2 M t Mf Mi 

JlJ Mf — Mf Mj 3 Mf - Mf Mj 

and where Mj are the neutrino mass eigenvalues. In [17, 19] it was shown that several of 
the most constraining observables, notably the Z decay to charged leptons and sin 2 6 [62], 
depended on the combination: 



(NN') ee (NN')^ - 2aT ~ 1 - |d e | 2 - \9^\ 2 - 2 aT. (40) 


Since from Table III |# e | 2 + \9^\ 2 ~ 10 3 , 2 aT must be of similar order so as to be competitive 
with the tree contribution. From Eq. (38) 


2 aT ~ 


aA 2 \9 a \ 4 
16vr’ 


(41) 
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Arrv[GeV] 

FIG. 9: T parameter versus 1-loop correction to m v for different values of the L -violating param¬ 
eters fjLi and /U 3 . 


where A is the mass scale of the heavy neutrinos and 6 a /y /2 their mixing with the flavour 
states from Eq. (9). Thus, in order for 2 aT ~ \6 a \ 2 it is necessary that A ~ 10 — 100 TeV. 
And, since \0 a \ 2 ~ |y Q | 2 u| w /2A 2 ~ 10 -3 , then \Y a \ ~ 1 — 10, on the very limit of perturba- 
tivity but, a priori, an interesting possibility. 

Furthermore, notice that the second term in Eq. (38) has the typical structure in the 
elements of the mixing matrix U of L-violating processes, such as, for example, neutrinoless 
double P decay. Indeed, this term stems from the correction to the Z propagator with two 
neutrinos running in the loop and a Majorana mass insertion and it is easy to see that it 
vanishes in the limit of exactly conserved Lepton number, taking all e* and fij to zero. Thus, 
if B — L is approximately conserved, the first term in Eq. (38) dominates the contribution 
to T. However, it can be shown that the matrix f,j is positive semi-definite for three extra 
heavy neutrinos or less 4 and can then be diagonalized as VAAfcU^, where V is a 

Unitary matrix and > 0. Thus, if B — L is approximately conserved: 

2 

> 0 . 

But from Eq. (40) T < 0 is mandatory so as to have the cancellation between T and \0 a \ 2 
discussed in [24], Thus, significant violations of B — L are necessary so that the second term 
in Eq. (38), which is allowed to be negative, can dominate over the first. 

Notice that, for arbitrary values of the B — L-violating parameters e* and /j,j, Eq. (8) is 
a completely general parametrization of a type-I Seesaw mechanism with three extra right- 
handed neutrinos. But, given Eq. (12), only jii and /i 3 are allowed to be sizable given the 



aT 


a 


E 


Sjrs w M w „ gj 




Preliminary explorations indicate that this argument can be generalized to more extra heavy neutrinos. 
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present constraints on the light neutrino masses and mixings. If \/ii\ A, A',/Z 3 a negative 
T can indeed be obtained: 


T ~ 


7 ew 






(43) 


If both n 1 and /13 are simultaneously included and dominate over the L-conserving A and 
A' then T is given by: 


T ~ 




4 

EW 


64v ts^M^ 



6/i/-h - (3^? + / 1 2 ) log 

^ 3 /ii 


(44) 


where /i = y//If + 4 / 13 . In this limit, negative values of T are also easily accessible. However, 
the price to pay is high, the approximate B — L symmetry protecting the Weinberg operator 
despite the Yukawas at the very border of perturbativity and the low Seesaw scale, has 
been strongly broken by /ii and /13. While this does not induce any dangerous corrections to 
neutrino masses at tree level, and hence when working with the Casas-Ibarra parametrization 
as in Ref. [24] the correct masses and mixings seem to be recovered, loop corrections need 
to also be taken into account since no protecting symmetry can now suppress them. Indeed, 
the loop contributions mediated by //] and /i 3 to the light neutrino masses are found to 
be [43, 63-66]: 

Am “"» = ( 3M ^ /(Mz) + > ( 46 ) 


with: 


(/i + /ii) 2 log(^r) _ (A-/^i) 2 iog (w) 

(/i + /ii ) 2 — 4 M 2 (/i — /ii ) 2 — AM 2 


(46) 


These corrections can indeed be sizable and in Fig. 9 we show the values that the loop 
contribution to the light neutrino masses take in order to recover a given value for —2 aT for 
different values of Hi and /r 3 . From inspection of Eq. (46), the limit of vanishing /ii would 
render f(M ) = 0, keeping under control the loop corrections to neutrino masses 5 . However, 
from Eq. (44), |/xi| > 1.91/x 3 1 is necessary for T < 0. Indeed, as can be seen in Fig. 9, if 
—2 aT ~ 10 -3 so as to implement the cancellation between tree and loop level contributions, 
corrections to the light neutrino masses ranging from ~ 100 keV to ~ 100 MeV, far exceeding 
present constraints, would be obtained. Thus, we conclude that, while the qualitatively 
important cancellations described in Ref. [24] can in principle take place and affect the 
constraints on the heavy-active neutrino mixing for Y a ~ 1 and A ~ 10 TeV, in practice 
large violations of the protecting B — L symmetry would be required, leading to too large 
radiative corrections to light neutrino masses. 


5 In this limit with ^3 A, A', L-symmetry is recovered with two degenerate neutrinos with mass /r 3 that 
form a Dirac pair. Hence, the symmetry ensures the stability of v masses at loop level but conversely 
drives T to positive values. 
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V. CONCLUSIONS 


In this work we have analyzed in detail the importance of loop corrections when deriving 
constraints on the mixing between the SM flavour eigenstates and the new heavy neutri¬ 
nos introduced in Seesaw mechanisms. Although naively the expectation is that radiative 
corrections involving these new states would be irrelevant given their weaker-than-weak 
interactions due to their singlet nature and, a priori, suppressed mixings with the SM neu¬ 
trinos, Seesaw models allow Yukawa couplings to be sizable, even order one. Thus, loop 
corrections involving Yukawa vertexes, when the loops involve the heavy neutrinos and the 
Higgs or the W and Z Goldstones, can indeed be sizable as shown in Ref. [24], In that work, 
it was shown that, for the low scale Seesaw mechanisms characterized by large Yukawas 
and low (electroweak) Seesaw scale, the contribution of the new degrees of freedom to the 
oblique parameters could indeed become as important as the tree level effects in some re¬ 
gions of the parameter space. Moreover, it was observed that several observables shared a 
common dependence between the T parameter and the tree level contribution, stemming 
from the modification by these effects of the muon decay through which Gp is determined 
and subsequently used as input for other observables. Thus, a partial cancellation between 
these tree and loop level contributions can significantly relax the bounds derived from these 
observables. Indeed, in Ref. [24] a good fit with sizable mixing was obtained in which the 
most stringent limits were avoided through this partial cancellation while standing tensions 
between the SM and some observables like the invisible width of the Z were alleviated. 

We have extended the analysis performed in Ref. [24] to include also vertex corrections 
and not only oblique parameters, since the sizable contributions from the heavy Yukawas 
do not vanish when taking the light neutrinos and charged lepton masses to zero. We find 
that, all in all, the oblique parameters do tend to dominate over the other loop corrections 
and their contribution could be sizable in some part of the parameter space. However, our 
MCMC scan shows no preference for any sizable loop corrections and the partial cancellation 
found in [24] is not reproduced. We then studied in detail the values of the T parameter 
preferred by data through our MCMC scan and saw that they were not only negligible, but 
always positive in our results. Indeed, for the cancellation between tree level contributions 
and the T parameter to take place, the latter must have negative values. We thus studied 
the necessary conditions for sizable negative values of the T parameter and realized that, not 
only sizable Yukawas and relatively low Seesaw scales are required, but also large violations 
of B — L. We then identified the only parameters in the mass matrix with three extra heavy 
neutrinos that could provide the necessary B — L violation required for T to be negative and 
competitive with tree level contributions, while keeping neutrino masses within their current 
bounds despite the large Yukawas, low Seesaw scale and loss of protecting B — L symmetry. 
Finally, we studied how these parameters would contribute to neutrino masses at loop level 
and found that, for the size of T required for the cancellation to take place, light neutrino 
masses would range from 10 keV to 100 MeV, effectively ruling out this possibility. 

We conclude that loop level corrections are only relevant in a small fraction of the Seesaw 
parameter space characterized by large Yukawa couplings and low Seesaw scale and that 
these corrections tend to strengthen the tree level contributions unless large deviations from 
B — L are present. If B — L is approximately conserved, data thus prefer regions of the 
parameter space where these loops are irrelevant. On the other hand, if B — L is strongly 
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violated, the cancellation discussed in Ref. [24] can indeed provide a good fit to data with 
a very relevant role of the loop contributions. However, these large violations of B — L at 


loop level also lead to too large contributions to the light neutrino masses and hence this 
possibility is ruled out. We therefore conclude that loop corrections can safely be neglected 


in analyses of the heavy neutrino mixings in Seesaw mechanisms. 

Finally we have also obtained relevant constraints on this mixing when B — L is an 
approximate symmetry, so as to recover the correct neutrino masses and mixings observed 
in neutrino oscillation searches. We find a mild (~ 90% CL) preference for non-zero mixing 
with the e flavours with a best fit at 9 e = 0.0341^014 or = 0.0351 o;q?4 for normal and 
inverted mass hierarchy respectively. In the case of normal hierarchy, this preference also 
induces non-zero mixing with the r flavour 6 r = 0.0181qq]J so as to recover the correct 
pattern of neutrino masses and mixings. On the other hand, small 6 ' M is preferred so as to 
keep n —y ey at acceptable levels in presence of non-zero 9 e . At the 2a level the following 
upper bounds are found: 9 e < 0.051, 9 M < 0.037 and 9 T < 0.049. 
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Appendix 


In this Appendix we list the self-energies, counterterms and diagrams that enter in the 
renormalization of the observables studied in Section III. 

Lepton-flavour-dependent counterterms: <5^ T n and <5 CT z 

The unrenormalized charged lepton fields Z° a can be written in terms of the renormalized 

lLa ones as 



( 47 ) 
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The most general expression for the Ip —> l a transition amplitude between fermionic Dirac 
states can be written as follows: 


M = (p 2 ) + fPR^f (p 2 ) + Pi^Sp (p 2 ) + P a E% (p 2 ) , (48) 

where E L = E L t and T, B = l n the physical observables only the Hermitian part of 5Z l 
appears and it is given by 



= - (rn 2 ) - mg mg(Y%g (mg) + (mf) ) + 


(sSK) + sf"K))], 

(49) 


with S'(p 2 ) = dE (p 2 ) /dp 2 . Therefore, the heavy neutrino contribution to SZ lcp can be 
obtained simply computing 

0 ± 


/ \ 





->— 

: 

P 
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MiUg t U’ ak 


(■ P L mp + P R m a ) B 0 (p 2 , Ml , M^) 



m a mp 

~Tif 



B\(p 2 , Ml, Mw) 


1 


(50) 


where Bi (and later Bij, Cij, Di and D t j) are the Passarino-Veltman integrals [67] using the 
notation from Ref. [68]. 

Similarly, the unrenormalized neutrino fields z/£- can also be written in terms of the 
renormalized ones as 

U Li — 0>ij + \ VLj- (51) 

The transition amplitude between two Majorana states reads 


E[f (f) = fP L Eg (p 2 ) + fP R Eg* (p 2 ) + P L E" (p 2 ) + P R E (p 2 ) , (52) 

where E L = E R * and E M = S M *. In the Majorana case, the Hermitian part of 8Z U can be 
written as 


iZT =5 + «?£) 


Sj (m?) - k (=« K 2 ) + S« ' (A) ) + R"' H) + E,f'(m 2 ) 


(53) 


Analogously to the charged lepton case, 5Z neu can thus be obtained from the heavy neutrino 
contribution to the neutrino self energy: 
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l - A V 
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167r,s^ v M^ fc=4 

5i(p 2 , M 2 , Mf) + P^p 2 , Ml Ml) 

+ |*Pr (MjCj k + M k C* k ) (MiCki + M k C* ki ) 
x 'bIpImImD + bIpImImD 
+ P L M fc (Mj-Qfc + M k C * k ) (M fc CW + MiCH) 


(54) 


x 


P 0 (p 2 , M 2 , Mf) - P 0 (P 2 , M 2 , M 2 ) 

+ PflMfc (Mjq k + M fc Q fc ) (M k C ki + MiC k i) 


x 


B 0 (pl Ml Ml) - B 0 (pl Ml Ml 


Finally, the lepton-flavour-dependent combinations that will correct and cancel the di¬ 
vergences of 1-loop corrections to the vertex Wul a and Zuu are respectively: 


U n 


e sz%uh + V ursz, 

i= 1 \/3=l J=1 

6 

■5° T 2 = E pzrc w + SZ^C lk ) 


neu 

ij 


(55) 

(56) 


k= 1 


Vertex interferences: V'l and V? 



h, 0°^-< 


E 3 nn*rji\Z 

i=l 1 o J-Wvila 


v w =- 

v rx - 


a 


Eti l T ol 2 

3 6 

EE M 2 k U ai U* ak C* ki [Coo(0,0, Ml Ml Ml) + C 00 ( 0,0, Mf, M 2 , M^)], 

( 57 ) 
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up to higher order corrections and where Tq is the corresponding tree level amplitude. 
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+ 2(C 00 (0, Ml, Ml, M fc 2 , M r 2 ) + C 00 (0, M 2 Zl M h 2 , Ml M r 2 ))" 


E [4C fei a fc Ml(C 00 (0, Ml, Ml, Ml, Ml) + C 00 (0, Ml, M fc 2 , Ml, Ml)) 


/c=4 


(58) 


up to higher order corrections and where T 0 is the corresponding tree level amplitude. 


Box contribution to p decay: B a p 
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+ 3(B 2 (M=) + D 2 (M|)) + 2 (D 3 (A© + D 3 (M|)) 


( 59 ) 


up to higher order corrections and where T 0 is the corresponding tree level amplitude and 
using the simplified notation Dij(M 2 ) —> Dij(0, 0, 0, M 2 , M 2 , M%, M^). Apart from the 
explicit sum over final state neutrinos in Eq. (59), the integral over the phase space is to be 
understood in both the numerator and denominator. 


Flavour-universal corrections to the gauge boson propagators: 


runiv N 
°W 


and 5 


univ N 

Z 


We label T,ww and Zzz the terms proportional to in the W and Z self-energies 
respectively. Notice that the SM contribution has been subtracted from the total self-energy, 
as we are interested in the contribution stemming from the new extra neutrinos. 
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x (2S 00 (p 2 , M 2 , M 2 ) + p 2 (^(p 2 , Ml Mf) + B u (p 2 , Mf, M 2 )) 


( 61 ) 


-3 


2-Boo (p 2 , 0, 0) + p 2 (-Bi(p 2 , 0, 0) + Bn(p 2 , 0, 0)) 


Notice that both in Eq. (60) and in Eq. (61) the sums run over all neutrino mass eigenstates 
(heavy and light) so here M* can represent both the heavy or the light neutrino masses. 

The oblique universal corrections to the electroweak observables can be written as a 
combination of the three following independent parameters [30, 31]: 
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and the renormalized self energies are given by: 
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( 66 ) 


Notice that, in the on-shell renormalization scheme Ef w (M'w) = E {Mf) = Ef (0) = 
Ef (0) = 0. Moreover, there is no contribution to the propagator of the photon from the 
extra heavy neutrinos and therefore Ef and Ef can be set to zero in the previous equations. 
In addition, there is no correction to E^ 7 either, so that Ef can be set to zero too. The 
universal oblique counterterms presented in Sec. Ill can thus be written as: 
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